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We provide a semiclassical theory of tunneling decay in a magnetic field and a three-dimensional potential 
of a general form. Because of broken time-reversal symmetry, the standard WKB technique has to be 
modified. The decay rate is found from the analysis of the set of the particle Hamiltonian trajectories 
in complex phase space and time. In a magnetic field, the tunneling particle comes out from the barrier 
with a finite velocity and behind the boundary of the classically allowed region. The exit location is 
obtained by matching the decaying and outgoing WKB waves at a caustic in complex configuration space. 
Different branches of the WKB wave function match on the switching surface in real space, where the 
slope of the wave function sharply changes. The theory is not limited to tunneling from potential wells 
which are parabolic near the minimum. For parabolic wells, we provide a bounce-type formulation in a 
magnetic field. The theory is applied to specific models which are relevant to tunneling from correlated 
two-dimensional electron systems in a magnetic field parallel to the electron layer. 
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I. INTRODUCTION 

Magnetic field can very strongly change the tunnel- 
ing rate of charged particles. This change, in turn, 
strongly depends on properties of the system, as in the 
well-known effect of giant hopping magnetoresistance in 
solidsEJ. Therefore tunneling in a magnetic field has been 
broadly used as a sensitive and revealing probe of election 
systems in solids, including quantum Hall systemsocl, 
two-lawes heterostructures away from the quantum Hall 
rcgionETu, and correlated electron layers on the surface 
of liquid hcliuni30. 

Correlated two-dimensional (2Dl-electron systems are 
currently attracting much interestlia. The possibility to 
extract information about electron correlations and dy- 
namics through tunneling in a magnetic field is one of the 
motivations of the present work. The major motivation, 
however, comes from the fact that tunneling in a mag- 
netic field is an interesting and in many respects unusual 
theoretical problem, even in the single-particle formula- 
tion. Existing results, although often highly nontrivial, 
are limited to, the cases wher.e-|tlie potentiaLias either a 
special formlijilj (e.g., linearis or parabolidi3) , or a part 
of the-notential or the magnetic field are in some sense 
weakEra 

The problem of tunneling has two parts. One is to find 
the tail of the wave function of the localized intrawell 
state under the potential barrier U(r) and behind it, and 
the other is to find the escape probability W. In the mag- 
netic field, W differs exponentially from the probability 
to reach the boundary of the classically allowed range 
U(r) = E, where E is the energy of the particle. This is 
because, as it tunnels, the particle is accelerated by the 
Lorentz force, and it comes out from the barrier with a 
finite velocity v. The standard argument that the exit 



point is the turning point v = relies on time-reversal 
symmetry (see below) and does not apply in the presence 
of a magnetic field. 

A simple potential U(r) and the wave function ip(r) 
of the metastable state in this potential are sketched in 
Figs. |l|, 0. The wave function decays away from the 
potential well. At some point r, on the background 
of the decaying tail there emerges a propagating small- 
amplitude wave packet, which corresponds to the escaped 
particle. As a result, in a part of the classically allowed 
region U(r) < E the function ip( r ) is determined by this 
wave packet, whereas in the other part of this region ip is 
determined by a different branch of the tail of the local- 
ized state. The boundary between these areas is sharp, 
and the slope of the wave function changes on this bound- 
ary nearly discontinuously. 

An important part of the WKB formulation of the-tim- 
neling escape problem in a magnetic field was foundEJ^ 
in the analysis of decay for a special model of an atomic 
system [see Eq. (||) below]. In a general case, both the 
tail of a metastable state and the outgoing wave packet 
can be-Jjound using the approach briefly outlined in our 
LetteiH 

In the WKB approximation the wave function is sought 
in the form 

ip(r) = D(t) exp[iS(r)] (fi = 1). (1) 

Here, S(r) is the classical action, and D is the prefac- 
tor. In the classically allowed range, ([[]) describes a wave 
propagating with a real momentum p = V5. On the 
other hand, in the classically forbidden range the wave 
function decays. For the ground intrawell state, the decay 
of ip( r ) is n °t accompanied by oscillations in the absence 
of a magnetic field. Then the action S(r) is purely imag- 
inary under the barrier, and | VS*! is the decrement of the 
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wave function. 

Both behind and inside the barrier, the action can be 
obtained from the Hamilton equations of motion 



-it, p — > ip, r 



U(r) -> -C/(r). 



(4) 



S 1 = p • r, r = dH/dp, p = —dH/dr, 
where if is the Hamiltonian of the system, 
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(A(r) is the vector potential of the magnetic field). 
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FIG. 1. Tunneling in a two-dimensional potential U(x,z) 
transverse to a magnetic field B pointing in the y direction. 
Initially the particle is localized in a metastable state behind 
the barrier (on the small- z side), with energy E. In contrast 
to the case B = 0, a particle emerges from the barrier with 
a finite velocity, and therefore the exit point is located away 
from the line U(r) — E. 




FIG. 2. The absolute value of the tunneling wave function 
(schematically). The maximum is located inside the potential 
well, i.e. behind the barrier in Fig. |l|. A Gaussian wave packet 
of nearly constant height describes the escaped particle, which 
shows up with a finite velocity. 

In the standard approach to tunneling decay, which ap- 
plies for B — OiJ'EacZl, one looks for the purely imaginary 
action S under the barrier. It is calculated by changing to 
imaginary time and momentum in Eqs. (^), but keeping 
coordinates real, 



Eqs.(||) then take the form of equations of classical mo- 
tion in an inverted potential —U(r), with energy — E > 
— U(r). The imaginary-time trajectory goes from the 
turning point inside the potential well to the turning 
point on the boundary of the classically allowed region, 
where it matches the appropriate classical trajectory of 
the escaped particle behind the barrier, with real r, p, t. 
The velocities on the both trajectories at the matching 
point can coincide only if v = 0. 

In the presence of a magnetic field, because of bro- 
ken time-reversal symmetry, the replacement (||) may not 
be performed. It would lead to a complex Hamiltonian, 
which makes no sense and indicates that a more general 
approach is required. The action S(r) is complex under 
the barrier for real r. This complexity has also an im- 
portant counterpart in the instanton formulation of the 
problem of tunneling decay in a magnetic field, see be- 
low. We note that a complex action arises also in-other 
cases, like barrier penetration for oblique incidences and 
scattering by a complex potential (as in the case of an 
absorbing medium)cj. The method discussed below can 
be applied to these problems as well. 

In this paper we consider a single-particle tunneling 
decay in a three-dimensional potential of a general form, 
for arbitrary magnetic fields. We illustrate the approach 
using a toy model of a correlated 2D electron system. 
We show that the tunneling exponent S(r) and the es- 
cape rate in a magnetic field can be found from dynamical 
equations (0) by analytically continuing these equations 
to complex phase space and time. The initial conditions 
for the trajectories are determined by the analytical con- 
tinuation of the usually known intrawell wave function. 
The resulting set of complex trajectories has singularities, 
caustics, on which there occurs branching of the complex 
action S(r) and the tails of the decaying and propagat- 
ing waves are matched. Careful analysis allows us to find 
the complete semiclassical wave function and reveal the 
singular features of ip( r ) related to the branching of S. 

In Sec. H and Appendix ^ we provide a simple model 
which catches basic physics of tunneling from correlated 
2D electron layers. In Sec. |l| we consider the tunnel- 
ing exponent and formulate the boundary value problem 
for tunneling trajectories in a magnetic field in complex 
phase space. In Sec. IV we discuss matching of differ- 
ent semiclassical solutions across the caustic of the set 
of the tunneling trajectories. We show that a switching 
surface (one of the anti-Stokes manifolds) starts at the 
caustic. The wave function has an observable singular 
feature at this surface, which is a sharp change of the 
slope of In |V>( r )l- ln Sec. [V| we provide explicit results 
for two simple exactly solvable models of physical inter- 
est, which also illustrate general features of tunneling in 
a magnetic field. In Sec. VI we discuss the path-integral 
formulation of t he p roblem of tunneling decay in a mag- 
netic field. Sec. VII contains concluding remarks. 



2 



II. A MODEL OF THE TUNNELING BARRIER 
FOR A CORRELATED 2D ELECTRON SYSTEM 

One of the most interesting systems where tunneling 
in a magnetic field has been investigated experimentally 
is a correlated 2D electron system. Here, electrons are 
localized in the z-direction in a metastable ID potential 
well Uq(z). The intrawcll electron motion is quantized in 
the z-direction, and electrons can tunnel from the well 
into extended states. Many 2D systems of current inter- 
est are strongly correlated: electrons are far away from 
each other, exchange is weak, and there is at least short- 
range order in the (x, y)-plane. The tunneling electron 
can be then identified and "labeled" . Its tunneling mo- 
tion is accompanied by the motion of other electrons. 
The many-electron dynamics of a correlated system can 
be described in terms of in-plane electron vibrations, and 
the corresponding Hamiltonian is given in Appendix |A] 
assuming that the electrons form a Wigner crystal. Here 
we will made a further simplification and think of an elec- 
tron as tunneling in a static potential created by all other 
electrons. As we showed earlier, this is a good approxi- 
mation for the tunneling problemliil. 

The static potential from surrounding electrons can 
be assumed to be parabolic with respect to the in-plane 
coordinates x, y. If the characteristic width L of the tun- 
neling barrier is less than the interelectron distance, the 
overall potential is a sum of the parabolic in-plane part 
and Uq(z), 

U(r) = ^(x 2 +y 2 ) + U (z). (5) 

The form of Uq(z) depends on the system. Inside the 
well Uo is often singular, like in the case of electrons on 
helium where Uq includes the image potential. The po- 
tential barrier itself can be close to a square barrier, as in 
the case of unbiased semiconductor heterostructures, or 
can be nearly linear, as in the presence of strong enough 
bias voltage, with 



U (z) = 




inside the barrier. The potential (g) describes, in par- 
ticular, the barrier for a correlated 2DES on a helium 
surface. This system was experimentally investigated in 
Ref.tj, and showed an unexpected dependence of the tun- 
neling rate, oh B and electron density that we recently 
addressedlilcj. 



III. THE TUNNELING EXPONENT 

For a smooth tunneling barrier U(r), the underbarrier 
wave function ip{r) (H) can be obtained from the tun- 
neling trajectories (BT The initial conditions for these 
trajectories are determined by the tail of the intrawell 



wave function. They can be obtained even if the poten- 
tial U(r) is singular within the well, as in the case of an 
image potential or a stepwise potential in a semiconduc- 
tor heterostructure. 

To obtain the initial conditions for the trajectories (g), 
we can take a surface £ close to the well and yet in the 
range where U (r) is already smooth. The wave function 
ip(j) on £ is presumably known from the solution of the 
Schrodinger equation inside the well and is semiclassical. 
Only the exponent of this wave function is needed to find 
the initial conditions for (||), which take the form 

r(0) = r| s , p(0) = -i[Vln^(r)] s , (7) 

with the action 5(0) = — i pn^(r)] E . Only the lowest- 
order terms in h should be kept in the expressions for 
p(0), 5(0). The final result should be independent of the 
choice of E. 

The trajectories (^) with the initial conditions (j^) form 
a two-parameter set, in the case of 3D tunneling. The two 
parameters are the initial coordinates on the surface £. 
We can choose curvilinear coordinates [x\, X2, X3) so that 
^3 Is == 0, and respectively 2:3(0) = 0. The trajectories 
are then parametrized by xi(0), 2:2(0). 

To illustrate these arguments we consider the initial 
conditions for an electron with the potential (||), 
which tunnels from a 2D layer. Inside the metastable 
potential well the electron motion separates into a quan- 
tized motion in the normal to the layer ^-direction and 
in-plane vibrations. To slightly simplify the analysis, we 
will neglect the effect of a magnetic field on the intrawell 
wave function, but not on the tail of ip(t) deep under 
the barrier where the effect will have accumulated. This 
is a good approximation for not too strong fields pro- 
vided the characteristic intrawell localization length 1 /*y 
is small compared to the tunneling length L. 

It is convenient then to choose the surface S as a plane 
z = const close to the well, but behind the intrawell 
turning point. We set z — 2:3 = on £ and choose 
x\ = x,X2 = y. If we set the energy of the out-of-plane 
motion E = and assume that the electron is in the 
ground intrawell state, we obtain from ([s|) 

z(0) = 0, Pz (0) = i% 5(0) = imLu [x 2 (0) + y 2 (0)]/2, 
p x (0) = imuj o x(0), p y (0) = imu Q y(0). (8) 

To find the trajectory which arrives at a given real r 
deep under the barrier, it may be necessary, particularly 
for B =/= 0, to start with a complex r(0). The correspond- 
ing values of p(0) can be found by analytically continuing 
the intrawell wave function to complex r. The whole tra- 
jectory will then go in complex phase space (including the 
configurational space) and also in complex time t. The 
energy of the trajectory is given by the energy of the in- 
trawell state from which the system tunnels. It remains 
real. 

The rate of tunneling decay is determined by Im 5 
at the point where the particle emerges from the barrier 
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as a semiclassical wave packet. This wave packet prop- 
agates along a real classical trajectory r c \(t), which is a 
real-time solution of the Hamiltonian equations (Q) . The 
underbarrier trajectory for tunneling escape should co- 
alesce with this classical trajectory. Therefore at some 
time t it should have both real coordinate and velocity, 



Im r(t) = Im p(t) = 0. 



(9) 



Eqs. (^|) determine the complex starting point of the 
trajectory for tunneling escape r(0) [i.e., the complex 
£1,2(0), since 2:3(0) = 0] and also the imaginary part 
of the duration of motion along this trajectory Im t. The 
real part of t remains undetermined: a change in Re t in 
(0) results just in a shift of the particle along the clas- 
sical trajectory r c \(t), see Fig. |[ Such a shift does not 
change Im S, since p = V5 is real along r c i(£). We note 
that the number of equations (J^j is equal to the number 
of variables Re £1,2(0), Im £1,2(0), and Im t [the value 
£3(0) = is fixed on £], with account taken of-energy 
conservation. The conditions (^|) were first giventS b ) for 
a (5-shape potential well and a linear tunneling barrier, 
but only the condition Im r(t) =0 was used. 
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FIG. 3. (a) Complex t plane for integrating the Hamilto- 
nian equations (^) in the escape problem. The line Im t = 
const corresponds to the classical trajectory of the escaped 
electron, which is shown in (b). Bold solid lines in (a) and 
(b) indicate the range where the amplitude of the propagat- 
ing wave exceeds the amplitude of the decaying underbarrier 
wave function. The escaped particle shows up as a semiclas- 
sical wave packet, with a finite velocity, at the point (full 
circle) where the classical trajectory intersects the switching 
line [thin solid line in (b)] . The crosses mark the value of t (a) 
and the position (b) of the caustic where it goes through real 
space. For the chosen parameter values, the time when the es- 
cape trajectory hits the caustic is numerically very close to the 
position of the cross in (a) . The specific data refer to tunnel- 
ing through the potential barrier (g), (^) transverse to a mag- 
netic field, which points in the y-direction, with UqTq = 1.2 
and lu c to — 1.2; time in (a) is in the units of to = 2mL/ / y. 

In the absence of a magnetic field, we can choose the 
surface £ such that the momentum p|s is imaginary for 
real r, i.e. the decay of the localized wave function is not 
accompanied by oscillations [cf. (||)]. The equations of 
motion (^|) can then be solved in purely imaginary time, 
with real r(t) and imaginary p(t). The escape trajectory 



ends at the turning point p = 0, even for a multidimen- 
sional systemEH. 

The tunneling exponent R is given by the value of Im S 
at any point on the trajectory r c i, 



R = 2 Im S*(r cl ). 



(10) 



For a physically meaningful solution, Im S should have a 
parabolic minimum at r c i as a function of the coordinates 
transverse to the trajectory. The outgoing wave packet 
will then be Gaussian near the maximum. 

From (^), even in the presence of a magnetic field the 
tunneling exponent can be obtained by solving the equa- 
tions of motion (|^) in imaginary time, with complex r. 
However, such solution does not give the wave function 
for real r between the well and the classical trajectory 
r c i. Neither does it tell us where the particle shows up 
on the classical trajectory. 



IV. BRANCHING OF THE ACTION AND ITS 
OBSERVABLE CONSEQUENCES 

The complete WKB solution of the tunneling problem 
can be obtained and the wave function ip(r) can be found 
if one takes into account that the action S as given by 
Eqs. (H) is a multivalued function of r, even though it is a 
single-valued function of t and £1,2(0). This means that 
several trajectories (||) with different t and £1,2(0) can 
go through one and the same point r. However, except 
for the points on the switching surface (see below), only 
one of the branches of the action S(r) contributes to the 
wave function ip(r). 



A. Caustics in a magnetic field 

In multidimensional systems, branching of the semi- 
classical action generally occurs |©»-eaustics, or envelopes 
of the Hamiltonian trajectoriest^til, see Fig. [|. Caus- 
tics are multidimensional counterparts of turning points 
familiar from the analysis of tunneling in ID systemsE3. 
The prefactor D(r) in the WKB wave function (0) di- 
verges at a caustic. In the case of ID semiclassical mo- 

— 1/2 

tion along the z-axis we have D oc p z , and D — > 00 
at the turning points Zt, which are given by the condi- 
tion p z = 0. The action is branching at turning points, 
S{z) — S{z t ) cx (z — zt) 3 / 2 . Its behavior near caustics in a 
multidimensional system is very similar (see below) , with 
z — z t corresponding to the distance from the caustic. 

Since neighboring Hamiltonian trajectories @, (0) 
touch each other on a caustic, the one-to-one correspon- 
dence between the coordinates £1, £2, £3 on the trajec- 
tory and the parameters t, £i(0), £2(0) breaks. Therefore 
the equation for a caustic has the form 



J(r) = 0, J(r) 



d(xi,x 2 ,x 3 ) 

d(£i(0),£ 2 (0),t)- 



(11) 
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The Jacobian J(r) can be related in a standard way 
to the prefactor -D(r), which in turn is determined by 
the first-order (in Ti) correction to the action — iS^\ 
D{r) = exp^Wfr)]. The equation for S^{r) can 
be obtained from the Schrodinger equation by seek- 
ing the wave function in the form ip = exp(iS'), with 
S = - iS (1 \ This gives 2vVSW = -Vv, where 
771V = VS(°' + (e/c)A. The vector v gives the velocity 
on the Hamiltonian trajectory (Q). Taking into account 
that vV5 (1) (r) = dS w /dt and that Vv = d1nJ(r)/dt, 
where the time derivatives are taken along the trajectory, 
we obtain 



D(r) = const x [J(r)] 



-1/2 



(12) 



It follows from Eqs. (11]), fllgj ) that the prefactor D di- 
verges on caustics, and the WKB approximation does not 
apply therecd. 
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FIG. 4. The set of tunneling trajectories (solid lines) with 
Re x = Im 2 = and the caustic (dashed line) for an effec- 
tively 2D tunneling problem. The data refer to the potential 
(0), (|^) and a magnetic field along the y-axis. The parameters 
are the same as in Fig. ^| The classical trajectory of the es- 
caped electron is the real-time continuation of the trajectory 
2 from the point Im x = shown by the open circle. The 



momentum p is real at this point, as explained in Sec. VA 
Other trajectories do not describe escape, since for them the 
necessary conditions Im i = Im x = are not satisfied. The 
caustic goes through real space at the point z — z c , Im x = 
(the intersection of the dotted lines). The trajectory of the 
escaped electron does not go through this point. 

There are both formal and physical distinctions be- 
tween caustics for tunneling trajectories with and with- 
out a magnetic field. For B ^ 0, the trajectories are 
complex, and Eq. ( |i"T| ) specifies a complex surface. This 
surface intersects the real 3D space along a line. In con- 
trast, for B — 0, because of time-reversal symmetry, tun- 
neling trajectories lie in real configuration_space. In this 
symmetric case Eq. ([ll]) specifies a surfaceEZI rather than 
a line in the 3D space. As we show below, this distinction 
leads to observable consequences. 



B. Local analysis near caustics 

The analysis of the wave function and branching of the 
action Sir) at complex caustics in the magnetic field is 
similar to that for caustics in real space, including turning 
points in the ID case23. Near a caustic, it is convenient 
to change to the variables x' , y' , and z' which are locally 
parallel and perpendicular to the caustic surface, respec- 
tively (we set z' = on the caustic). Since a caustic is an 
envelope of the Hamiltonian trajectories (|^), the normal 
to the caustic component of the velocity is v z / =0 for 
z' = 0. However, for B ^ the normal component of 
the momentum is not equal to zero. Therefore the wave 
function near a point r caust on the caustic can be sought 
in the form 



■0( r c 



r = e 



i *"caust ) 5 



(13) 



where p ca ust is the momentum at the point r caust (note 
that r caust is a 2D complex vector, z' caust = 0). We as- 
sume that the dependence of p ca ust on r caust (i.e., along 
the caustic) is smooth, and the dependence of <j> on r caust 
is much more smooth than on z'. Generally, p caU st is 
complex even where the caustic goes through real space, 
Im r caust = 0. Therefore the classical trajectory of the 
escaped particle does not go through the caustic, in con- 
trast to the case of zero magnetic field, cf. Fig. || 

Because v z > = 0, the equation for <j>(z'; r caust ), which 
follows from the 3D Schrodinger equation with a mag- 
netic field, coincides with the ID Schrodinger equation 
near a turning point 



1 d 



2m dz 



72 + U' z ,(r ca _ us t)z' 



{z';r caust ) = (14) 



[here, U', = dV /dz']. The boundary conditions to this 
equation are discussed below. 

The function <\> is singled- valued_Jt is given by a linear 
combination of the Airy functionsE3. For comparatively 
large \z'\ (but still close to the caustic) it becomes a linear 
combination of the functions 



= (*') 



-1/4 



exp 



3/2' 



(15) 



with a = a(r caust ) = (2/3)[-2mU' z ,} 1/2 . To make the 
functions uniquely defined, we have to make a cut 
on the complex z'-plane. Our choice of the cut is shown 
in Fig. ||. 

With account taken of Eq. (p"3|), we find that the action 
near the caustic is 

S{x', y\ z') « S(x', y', 0) + ( Pcaust W + az' 3/2 , (16) 

with an appropriately chosen branch of z' 3 ^ 2 . 

Another way to understand Eq. (jl^) is based on the 
analysis of the set of the Hamiltonian trajectories (|^). 
Because the caustic is an envelope of the trajectories 
and v z > = on the caustic, z' is quadratic in the incre- 
ments fei^O), St of the parameters of the set. Therefore 
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<foi,2(0),<$£ are nonanalytic in z', as is also the action S. 
Taking into account cubic terms in Sxi^{0), St we obtain 
@. The coefficients in S can be expressed in terms of 
the derivatives of S, r calculated along the trajectories 
over £1,2(0), t. 



C. Choosing the action branches 

Eq. ( |l4| ) describes how the WKB solutions, which cor- 
respond to different branches of the action S, connect on 
the caustic. Of interest for the problem of tunneling es- 
cape is the caustic where there are connected the tails of 
the intrawell wave function and the outgoing wave packet 
for the escaped particle. From (|lT|), on this caustic there 
is a point through which there goes the Hamiltonian tra- 
jectory (||), (Q) for escape, with the initial coordinates 
£1,2(0) given by the condition (||) of arriving, ultimately, 
at the classical outgoing trajectory r c i, cf. Fig. 0. 
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FIG. 5. The function <p(z') (|14|) for large \z'\ on the com- 
plex z'-plane perpendicular to the caustic surface. The dashed 
lines, argz' = (2ra+l)7r/3, show the Stokes lines where the ra- 
tio of the functions \w2/wi \ ( |l5| ) reaches its maximum or min- 
imum. The anti-Stokes lines (solid lines), argz' = 2nn/3, are 
the lines where \w\\ = |i02|. At these lines the ratio |w2/wi| 
changes from exponentially large to exponentially small with 
varying argz'. The coefficients are found from the radiation 
boundary condition following the Stokes prescriptions. The 
dotted line shows the branch cut. 

In general, a caustic can be thought of as a "mirror", 
which partly reflects the wave packet. The boundary con- 
dition to Eq. ( |l4| ) for the tunneling escape problem is the 
"radiation condition" . Very far from the potential well 
the solution of the full Schrodinger equation is a semi- 
classical wave packet moving in real space away from the 
well. Respectively, the wave function far from the well 
has a form ip( r ) 00 exp[iS(r)]. This solution has to be 
continued to the caustic, which means that there is a 



range of directions in the complex z'-plane not too close 
to the caustic (|az' 3 ^ 2 | 3> h) where the wave function 
is described by only one exponential exp[i5'(r)]. Away 
from this range, the wave function is a combination of 
two waves. This is again similar to a ID problem, where 
behind the turning point zt the wave function for real z 
is a propagating wave, whereas before zt there is a wave 
with an exponentially decaying amplitude coming from 
the intrawell state and the wave reflected back to this 
state. 

To connect the tails of the wave functions near the 
caustic in our problem one can use the explicit solution 
of Eq. ( [Til ) jsath the radiation boundary condition, as in 
the ID caseEl 

An alternative way to sec how the boundary condition 
works, which also allows us to reveal the specific feature 
of tunneling in a magnetic field, is to follow the trans- 
formation of the wave function (f> for |cvz' 3 ^ 2 | > li as 
the argument of z' varies by 2ir. This analysis 
on the notion of the Stokes and anti-Stokes lineal 
We count argz' off from arg a 2 / 3 . Then, for the choice 
of the cut in Fig. |^ and with the functions 101,2 given by 
Eq. (|l5|), the Stokes lines are the rays arg z' = (2n+l)n/3 
with n = 0,1,2. On the Stokes lines Rez' 3 ^ 2 = 
and the ratio |u>2/wi| is extremal (maximal or minimal). 
The anti-Stokes lines are the rays argz' = 2nn/3 with 
n = 0, 1, 2, where |u>2/wi| = 1. 

From the radiation boundary condition, there is a 
range A of argz' where <fi(z') is given by only one of 
the functions wi^fz'), not by a superposition of w\ and 
W2- This condition is physically meaningful provided 
the corresponding Wi is exponentially small compared to 
i/>3_i in a part of the range A. [Otherwise the condition 
4>[z') = const x Wi is not a limitation, in the WKB ap- 
proximation.] It follows from the analysis below that A 
includes one of the rays argz' = or tt, along which c/>(z') 
is oscillating as exp(iaz' 3 ^ 2 ) or exp(— iaz' 3 ^ 2 ). 

For concreteness, we will assume that A contains the 
semi-axis argz' = 0, and that between argz' = and the 
cut in Fig. I (Hz') = Cw 2 {z') oc exp[i|a|z' L where C 
is a constant (it is not incorporated into Fig. |5|). 

Since the function <fi(z') is single- valued, if we cross the 
cut in Fig. |^ by incrementing argz', 4>(z') becomes equal 
to —iCwi(z'). It remains exponentially small as argz' 
grows up to 27r/3, including the Stokes line argz' = tt/3. 
Then behind the anti-Stokes line at argz' = 27r/3, the 
function w\ becomes exponentially big. It is important 
that, on the Stokes line argz' = tt, one has to take into 
account the admixture to <f> of an exponentially small 
term oc W2(z'). This can be seen from the explicit solu- 
tion for (j). The need to incorporate this term can also be 
understood by noticing that, when we increment argz' 
by 2n, we have to recover the original asymptotic form 
of ip. This latter argument explains the coefficient at W2 
in Fig. ||. 

On the anti-Stokes line arg z' = 47r/3, the values of \wi \ 
and \w 2 \ become equal to each other, and <fi(z') is primar- 
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ily determined by W2 for larger arg z' . After arg z' crosses 
the Stokes line 57r/3, the exponentially small term wij-ki 
4>{z') disappears, according to the Stokes prescription^! 

The wave function which connects to the outgoing 
wave on the caustic is that of the metastable state. Since 
the asymptotic behavior of <j){z') for large \z'\ is given by 
W2{z') 1 the asymptotic behavior of the wave function of 
the metastable state within the range 27r/3 < argz' < 
Att/3 is given by wi(z'). We note that, for the radia- 
tion boundary condition, the switching between the wave 
functions occurs only on one of the anti-Stokes lines. 

D. Switching between the branches of the wave 
function 

Switching between the WKB wave functions of the lo- 
calized state and the outgoing wave packet is an impor- 
tant observable consequence of the analysis in the previ- 
ous subsection. It is due to branching of the WKB action. 
The switching manifold starts on the caustic and is given 
by the condition Im Si(r) = ImS^r), where Si 2 are the 
actions for the corresponding WKB branches. On the 
opposite sides of the switching manifold one of the WKB 
wave functions is exponentially bigger than the other. 

In the presence of a magnetic field, caustics go through 
real space along the lines given by the condition J(r) = 0, 
Im r = 0, with r = r(xi.2(0), i) being a point on the 
Hamiltonian trajectory (g), (|7|). The switching manifold 
in real space is a surface which starts from the caustic 
line and goes away from it in one direction. Although the 
wave function is continuous on this surface, the derivative 
of its logarithm sharply changes from VSi to VS2. 

The exit point r in the configuration space where the 
escaped particle emerges from under the barrier is deter- 
mined by the intersection of the classical escape trajec- 
tory r c i(i) and the switching surface. This point can be 
found from the global analysis of the WKB wave func- 
tion. It does not lie on the caustic, nor is it given by 
the condition p = or equivalently, U(r) = E. In a 
2D system the caustic pierces real space at a point, and 
the switching surface becomes a line. An example of the 
caustic, the switching line, and the exit point for a 2D 
system is shown in Fig. ^. 

In the absence of a magnetic field, as we mentioned 
before, the caustic for the tunneling trajectories is a sur- 
face rather than a line in real configuration space. The 
switching surface in real space should then coincide with 
the caustic. The exit point is the turning point p = 0. 
The whole topology thus differs qualitatively from that 
for 

V. TUNNELING FROM A CORRELATED 2D 
ELECTRON SYSTEM 

We now apply the general approach to tunneling from 
a correlated 2D electron system transverse to a magnetic 



field and illustrate the occurrence of the singularities dis- 
cussed above. We will use the simple but nontrivial 
model of the electron system discussed in Sec. |0| and 
described by Eqs. (||), and its generalization to the 
case where the in-plane symmetry is broken. We will as- 
sume that the magnetic field is parallel to the electron 
layer, and choose the j/-axes along the field B. 

A. A model with in-plane symmetry 

For an electron in the potential (||), (||), classical mo- 
tion along the B || y axis is decoupled from the motion 
in the (x, z)-plane. The WKB tunneling problem then 
becomes two-dimensional, with complex classical trajec- 
tories (|^) lying in this plane. The Hamiltonian equa- 
tions are linear, and we can find the trajectories explic- 
itly. We can_also explicitly find the tunneling exponent 
and analyzes 3 ] its dependence on the two dimensionless 
parameters lvqTq and w c to, where tq is the tunneling time 
in the absence of the magnetic field, tq = 2mL/j. The 
expression for the tunneling exponent also follows as a 
limiting case from the result of the next subsection. Here 
we will discuss the structure of the WKB action. 

The symmetry of the potential U{x,y,z) = 
U(±x, ±y, z) in (g) gives rise to a specific symmetry of 
the set of the Hamiltonian trajectories, 

t-y-t*, x^-x*, z^z*, S^>-S*, (17) 

and of the singularities of this set. In particular, the caus- 
tic where there are connected the outgoing wave packet 
and the intrawell wave function goes through real plane 
(x, z) at a point x caus t = on the symmetry axis. The 
z-coordinate of this point depends on the magnetic field, 
with z caust = L for B = 0. The form of the action for 
z < z caust near the caustic in the symmetry plane x — 
is shown in Fig. @. As seen from the inset, the slope of the 
action is dim S/dz > at z caust , in contrast to the ID 
case where the slope is equal to zero at the turning point 
z t . We note that the branches 1 and 2 are formed by 
the trajectories that go through real space at times -Im t 
being, respectively, smaller and larger than -Im t caU st for 
the same trajectories [i C aust is given by Eq. 

Using the explicit form of the trajectories (0) with the 
initial conditions (^), one can find the complex caustic 
z'(x, z) = near x = 0, z = z caust . It has the form 
z - z caust = iC'x with real C'/B > 0, cf. Fig. f| It is 
seen from Fig. ^ that the singular parts of Im Si, 2 be- 
have as T^caust — z) 3 / 2 near the caustic. Therefore we 
can choose the coordinate z' that gives the distance from 
the caustic as z — z caust — iC'x. The branching behavior 
near the caustic is then described by Fig. [|. The range 
7T < argz' < 2tt corresponds to real z and real positive 
x. 

Close to the caustic, only one branch of the action 
describes the wave function for negative x and real z 
(the upper half of the complex- 2/ plane). For positive 
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x, we should keep both branches, and depending on 
x, z, the WKB wave function is given by the branch 
with the smaller Im S. Near the switching line where 
Im Si = Im 62 the total wave function is given by a 
linear combination of the two WKB solutions. 
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1.508 
























J 2 
p- 


1.507 




2/7 
/ 1 




2 
















Im 


1.506 




.1.1. 










1.05 


Z n 1-08 ^ 





















D 




0.5 






1 Z 



z/L caus 

FIG. 6. Two branches of the action on the symmetry axis 
x — as a function of the tunneling coordinate z before the 
branching point, for the symmetric model and same pa- 
rameter values as in Fig. ^. The vicinity of the cusp Zcaust 
is zoomed in the inset to show that the upper branch is non- 
monotonic. Its extremum at z m lies on the classical trajectory 
of the escaping particle shown in Fig. ^(b). However, the par- 
ticle emerges from the barrier for z > z m and x 7^ 0. 

In Fig. |^ for z < z C aust, the action branch 1 describes 
the tail of the intrawell wave function, and the branch 
2 corresponds to the wave "reflected" from the caustic. 
The branch 2 is nonmonotonic in z for x = 0, with a 
minimum at z m < z C aust- By symmetry the mo- 

mentum component p x is real for x — 0, whereas p z — 
at z m . Therefore the point z = z m , x = belongs to the 
classical escape trajectory shown in Fig. ||(b). Again by 
symmetry, this is also the point where the escape trajec- 
tory comes closest to the well at z = 0. However, this is 
not the exit point for the tunneling particle in the con- 
figuration space. Indeed, the wave function at this point 
is determined by the branch 1, because Im Si < Im S2. 

Cross-sections of the action surfaces by planes 
z = const are shown in Fig. ^. For z < z caust , both 
branches Im Si ,2 are symmetrical in x. However, the 
branch 2 is nonmonotonic in x for z > z m . It has a local 
maximum at x = and two symmetrical minima. These 
minima lie on the classical trajectory shown in Fig. §(b). 
For z — z m , the maximum and the minima merge to- 
gether. We note that at this point Im S2 oc x 4 near the 
minimum. 

Behind the caustic in real space, z > z caus t, the action 
S{x, z) on one of the two branches is equal to — S*(— x, z) 
on the other branch (cf. Fig. Q). At their minima with 
respect to x, the values of Im S(x, z) are independent of 
z. These minima lie on the escape classical trajectory, as 
seen from the comparison with Im S for z < z C aust- The 
branches of the action describe the wave packets incident 
on the barrier from large z and the emitted wave packet. 
Only the latter is physically meaningful for the problem 
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FIG. 7. Cross-sections of the function Im S by the plane 
z = const in the case of tunneling in a symmetric potential 

(|), for (a) Z m < Z < Zcaust, (b) Z = Zcaust, and (c) Z > Zcaust- 

The parameter values are the same as in Figs. ^, []. The solid 
lines show the branches of Im S that determine the exponent 
of the WKB wave function. The minima of the branch 2 lie 
on the classical trajectory shown in Fig. ^(b) . 

As discussed above, switching between the action 
branches occurs for positive x where the branches of Im S 
cross each other. The resulting action, which determines 
the WKB wave function, is shown in Fig. Q by solid lines. 
The switching line thus obtained coincides near the caus- 
tic with the anti-Stokes switching line discussed in the 
previous section. 

The escaped particle can be observed as a semiclassical 
object if Im S2(r c i) < Im Si(r c i). It "shows up" at the 
point where the classical escape trajectory intersects the 
switching line, cf. Fig. |3|. The exit point is located for 
x > even though the potential (|^) is symmetric. This is 
a consequence of the symmetry-breaking by a magnetic 
field. 



B. A non-symmetric model 

The problem considered in the previous section for 
the quadratic in x, y potential U(r) (^) can be solved 
differently. The trick iai-3 to make a canonical trans- 
formation to the new coordinate p x and the conju- 
gate momentum —x. The kinetic energy then becomes 
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[m 2 u!QX 2 +p 2 ]/2m and is independent of the new coordi- 
nates and the magnetic field. The time-reversal symme- 
try is thus "restored" , and the problem is mapped onto 
the standard problem of tunneling in the 2D potential 
Uq(z) + (p x + muj c z) 2 /2m. 

The general method discussed in this paper is not lim- 
ited to potentials with these special properties. In this 
subsection we illustrate how this method works where 
variables do not separate. To this end, we consider a 2D 
problem of tunneling transverse to the field B || y in the 
potential 

U(x, z) = \-mujlx 2 + fixz + X ( 1 - ( z > °)> ( 18 ) 
2 2m V LJ 

which differs from the potential discussed earlier by the 
term \xxz. In the problem of tunneling from a 2D electron 
system, this term mimics the dependence of the tunneling 
potential on the displacements of neighboring electrons, 
see Appendix [a|. 
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FIG. 8. Cross-sections of Im S(x, z) by the plane z = const 
in the case of tunneling in the asymmetric 2D potential 
U(x,z) (|l|), for (a) z m < z < z C aust, (b) z = Zcaust, and 
(c) z > z ca u S t. The values of lo c to and luqtq are the same 
as in Figs. ^, ^, and the dimensionless asymmetry parameter 
4/xmL 2 /7 2 = 1/2. The solid lines show the branches of Im 5* 
that determine the exponent of the WKB wave function. The 
minima of the branch 2 lie on the classical escape trajectory 
shown in (d). The cross in (d) marks the branching point of 
the function Im 5* where the caustic goes through real space. 
The switching line (thin solid line) starts at the branching 
point. 

The term /ixz breaks the symmetry of the Hamilto- 
nian trajectories ([l7|). However, the Hamiltonian equa- 
tions (Q) are still linear, and we explicitly solved them. 
The results for Im S and the classical escape trajectory 



obtained using the initial conditions (^J) are shown in 
Figs. ||, ^|. Because of broken symmetry, the branching 
point of the action in real space, where the caustic of the 
set of trajectories @ goes through real plane (x,z), lies 
at ^caust 7^ 0. It is marked by the cross on Fig. ||(d). 
Its position depends on fj, and other parameters of the 
system. Similarly, the time where the caustic crosses the 
real space has both real and imaginary part, in contrast 
to the case /i = where it was purely imaginary, see 
Fig. | 

For /j ^ 0, the surfaces Im S(x, z) become asymmetric. 
The general structure of the solution, however, remains 
the same as in the case fi = 0. This can be seen by 
comparing the cross-sections of the action in Figs. ^ and 
^. In both figures, the cross-sections in (a), (b), and (c) 
refer to the planes z < z ciiUSt , z = z caust , and z > z caust , 
respectively. As in the symmetrical case, the branches 1 
and 2 for fi ^ are formed by the trajectories in complex 
time with -Im t being, respectively, smaller and bigger 
than -Im i caust . At the point a; caust , z ca , ust , the branches 
1 and 2 touch each other. 
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FIG. 9. Escape trajectories for the dimensionless asymme- 
try parameter 4fimL 2 /"f 2 — 0,0.5,1.5. The filled circles are 
exit points, and the open circles are the points where the 
trajectory is most close to the localized state. 



For z < z caust there occurs switching between the 
branches. It can be analyzed in the same way as for 
/i = 0. The WKB wave function is determined by the 
branch in Fig. || shown with the solid line. The switch- 
ing line starts from the branching point x caU st, Zcaust and 
goes in the direction of positive x. 

The branch 2 has two minima as a function of x for 
z m < z < Zcaust, which are asymmetric for fi ^ 0. How- 
ever, their depths, i.e. the minimal values of Im S on 
this branch, remain equal to each other and are the same 
in all cross-sections z = const. These minima lie on the 
classical trajectory along which the electron escapes. At 
z = z m they merge together, and Im S becomes quartic in 
x near the minimum. The value z m shows how close the 
escape trajectory comes to the localized intrawell state. 
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The classical trajectory becomes observable in configu- 
ration space once it crosses the switching line. The shape 
of the trajectory and the exit point for several values of 
the asymmetry parameter [i are shown in Fig. ^. The out- 
going wave packet is Gaussian near the maximum (Im S 
is parabolic near the corresponding minimum). 

By solving the Hamiltonian equations for the po- 
tential (|l8|), we obtained the following expression for the 



Here, p, = /itq /to is the dimensionless asymmetry 
parameter. The dimensionless frequencies vq = loqTq, 
v c = uj c tq, and v — + v 2 ) 1 ! 2 characterize the motion 
under the barrier, v 2 = ±v 2 /2 + */ v A /4 + p? . 




FIG. 10. The tunneling exponent J? as a function of the 
magnetic field and the asymmetry parameter /i in the model 
for luqto — 1.2. The function R is even in fj,. 

The tunneling exponent R depends on the interrelation 
between the in-plane electron dynamics, which is char- 
acterized by the frequency ujq, the cyclotron frequency 
uj c , the tunneling time To, and the asymmetry-induced 
mixing of in-plane and out-of-plane motions p,. The de- 
pendence of R on the asymmetry parameter and on the 
magnetic field is shown in Fig. |l0|. For oj c = p, = 0, 
we have T r d = 1, and R = AjL/S. The exponent R in- 
creases with the magnetic, field. A qualitative result of 
the in-plane confinements is that it eliminates the local- 
izationpttwfche linear potential Uq(z) due to the magnetic 
barrierEffiiil, i.e. the divergence of R for lo c tq > 1. The 
effect of the magnetic field on R becomes small for strong 
confinement, ujqTo 3> 1 and u>q 3> oj c . 



tunneling exponent (|l0|), 

R = 2 7 L[r rd + u K(T ld )] (19) 

Here, T r d is the imaginary part of the time to reach the 
classical escape trajectory [see Fig. ||(a)] in the units 
of that same time for B = /j, = 0, which is given by 
To = 2mL/"f. Along with the function k in (|l9|), r rc j can 
be found from the equation 



(20) 

I 

The asymmetry results in lowering of the tunneling 
barrier for B = and the corresponding increase of the 
tunneling rate. This can be qualitatively understood, 
since a displacement in the cc-direction with [ix < in- 
creases the effective force in the z-direction which pulls 
the electron away from the layer. For small asymmetry, 
the correction to R is quadratic in p. In the limit 
of a thin and high tunneling barrier for x — or soft 
in-plane vibrations, where o->o T o <C h tunneling is most 
likely to occur in the barrier, which is adiabatically pre- 
pared by the in-plane displacement x (the "completely 
adjusted" barrierEj). The expression for R takes a form 

i?=2 7 L(3VA 2 ) 1/3 (^or «l), (21) 

it depends on the in-plane frequency ujq and on /i nonan- 
alytically. The role of the asymmetry increases with the 
magnetic field, as seen from Fig. p^, 

In terms of comparison with the currently available 
experimental data on tunneling from a correlated many- 
electron systemEHl, of utmost interest is the situation 
where the asymmetry is small. The observed dependence 
of R on B did not show the divergence expected for an 
unconfined electron. The simple model LH) provides a 
qualitative explanation of the experimented. An excel- 
lent quantitative agreement, without adjustable parame- 
ters, was achieved by incorporating the cur-stature of the 
potential Uq(z) due to electron correlations!!!!. 

VI. THE PATH-INTEGRAL FORMULATION IN 
A MAGNETIC FIELD 

In the absence of a magnetic field, the problem of 
tunncling|-decay is often considered using the instanton 
tcchniquecjO. This technique applies if the potential 
well is parabolic near the minimum and thermalization 
inside the well occurs much faster than escape from the 
well (in the case of 2D electron systems, both conditions 



, , v q ( r rd cos v_ T rd - v_ 1 sin v_ r rd ) + v 2 cos v_ r rd - v_ (1 - fa + t'oT'rd) sin V_ T r d 

K{T ld ) = — : r 

pr {Vq COS V T r d — V Sin V_ T r d) 

_ (^pTrd ~ V+) COshiy + Trd + (1 ~ VQ + ^OTrd - Vp / 1 V+) shlh V + T ld 

/t 2 (^o cosh f-f.Trd + v+ sinh is+T ld ) 
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are often violatedEl) . Because the Schrodinger equation 
for metastable states has to be solved with the radiation 
boundary condition, the energies of these state acquire 
small imaginary parts, and so does the partition func- 
tion Z. The escape rate W for finite temperatures is 
simply related to Im Z, 



W w 2TIm Z/Re Z 



(22) 



(we have set fcs = 1). 

The partition function is given by-ilie integral over pe- 
riodic paths r(r) in imaginary time~3, 



Z 



(0)=r(/3) 



Vv(T)exp[-S E [r(T) 



(23) 



where /3 — T _1 and Se is the Euclidean action (the ac- 
tion in imaginary time) . It is real for B = and for real 
trajectories r(r). 

The general expressions ( p2] ) , (^3|) should also apply in 
the presence of a magnetic field. However, the Euclidean 
action for an electron 



Se = 



,0 


m 




/ dr 






Jo 


~2 





e , , dr 
[/ r -H-A r — 

c dr 



(24) 



is now complex. Therefore the standard of eval- 

uating the escape rate has to be revised, excep t for special 
symmetric cases like the one discussed in Sec. |V A| , where 
nnp ran change to new varibales in which Se becomes 
realllJll3. Although we will often refer to the action func- 
tional of the form (p4]), much of the results below apply 
also to a more general retarded Euclidean action, which 
is of interest for systems coupled to a bath. 

In the spirit of the WKB approximation, the path 
integral (|2^) will be evaluated by the steepest descent 
method. The extremal paths r(r) of the functional Se 
satisfy the equation 



m—j = VU(r)+t- 
dr z c 



dr 

d^ XB 



(25) 



The equation of motion (|25J) has to be solved with the 
periodic boundary condition r(0) = r(/3). Note that the 
sign of the potential has been inverted compared to the 
case of classical motion in real time. 

For low temperatures and for the potential U(r), which 
is parabolic near its intrawell minimum r wc n , Eq. (|2^) has 
a solution r(r) = r wc n, with Se = 0. It gives the real 
part of the partition function, see below. As in the case 
B = 0, the imaginary part of Z is determined by an- 
other solution of (p5f), which is of the bounce type. This 
solution, rb(r), starts near r we n, slides downhill in the 
inverted potential —U(r), and in time (3 comes back. For 
B = the corresponding path is a symmetric real tra- 
jectory, i"b(r) — r b (/3 — t), which bounces off the turning 
point r b (/3/2) = 0. 

For B ^ 0, because of broken time-reversal symmetry, 
the path r^(r) is complex, and the velocity along this 



path does not become equal to zero. The path is not 
symmetrical in time, because the replacement r — > — r 
changes Eq. (p5|). However, if we simultaneously change 
i — > —i, the equation remains unchanged. Therefore the 
bounce-type path has the symmetry 



rb(r) 



(26) 



An immediate and very important consequence of 
Eq. (p6| ) is that the value of Ss(rb) for the bounce- type 
path is real. This value gives the tunneling exponent, see 
below. 



A. The eigenvalue problem 

The prefactor in Z can be found by integrating over 
the tubes of paths around the extremal paths. It can be 
done by expanding Se in deviations from the extremal 
paths to the second order, and then expanding r(r)— r wo n 
and r(r) — r b (r) in the eigenfunctions ip n (r) of the ap- 
propriate eigenvalue problem, 



dT / F(r ) r')V„(T')=A„V»(r), 
4-(r,r') = 5 2 S E /6r i (T)dr j (T'). (27) 



Here, the derivatives of the action are calculated on the 
corresponding extremal trajectory r we n and rb(r), and 
periodic boundary conditions are assumed. The op- 
erator F is simplified for a non-retarded action (E), 
F 1] {t 1 t') = 5{t-t')%{t). 

For B = 0, the operator j- is Hermitian, with fjj = 
m,5ij(d 2 / dr 2 ) + d 2 U/dridrj, if the action is given by 
Therefore the functions t/>„(t) form complete and 
orthogonal sets for each extreme trajectory (p5|), and the 
eigenvalues A„ are real. 

For B ^ 0, the operator J- becomes non-Hermitian. 
For example, in the case of a uniform magnetic field 
in ffci(r) acquires an extra term i{e / c)tkij B j{d / dr) 
(ekij is the Levi-Civita symbol). Therefore some of the 
eigenvalues A„ become complex. The eigenvectors ij) n 
with different n are orthogonal not to each other, but to 
the eigenvectors <f> n of the Hermitian conjugate operator, 



f dr' Ft (t,t>„(t') 
Jo 



a;0„(t) 



Taking into account the symmetry ( p6] ) of the extremal 
trajectories, we find that 



F^t,t') = F(/3-t,P-t') = F*(t,t'). 



(28) 



The energy spectra for several complex Hamiltonians 
with the symmetry, which is similar to ( p8| ) (and was 
called the PT-symmetry) , were investigated|-earlier nu- 
merically and using the WKB approximatior£3. 
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The symmetry (|2§|) has several consequences. First, it 
shows that ip n (T) = a„0*(r), where a n is a constant. 
This means that, with proper normalization, the orthog- 
onality relation becomes 







dTip m (T)ip n (T) = 5 V 



(29) 



(here, we assumed that the eigenvalues are nondegener- 
ate; for degenerate eigenvalues, the condition (|9|) can be 
satisfied by choosing appropriate linear combinations of 
the eigenfunctions with same A„). 

It also follows from Eq. ([28| ) that, if ip n ( T ) i s an eigen- 
function of (|2?|) with an eigenvalue A„, then ip* n (/3 — r) is 
also an eigenfunction of the same boundary value prob- 
lem, but with the eigenvalue A* . This means that a part 
of the eigenvalues A ra in (27) are real, whereas another 
part are pairs of complex conjugate numbers. 

Pairs of complex conjugate eigenvalues are formed 
in the following way. For B = 0, all eigenvalues are 
real. With increasing B some eigenvalues approach each 
other, pairwise, while still remaining real. Eventually 
they merge, and for larger B become complex conju- 
gate, as described in Appendix |^. For ID Schrodinger- 
type equation with complex Hamiltonians such behavior 
with the vacjring control parameter was indeed observed 
numcricallyO. 



2. Eigenvalues for the bounce trajectory 



A specific feature of the eigenvalue problem (|27j) for 
the bounce trajectory rf,(r) at low temperatures is that 
one of the eigenvalues is Ai = 0. It corresponds to the 
eigenfunction i/' 1 (t) oc rb(i~). For B = 0, the vector 
function ipi(r) °c *b( T ) nas one zero f° r au components. 
Therefore, as can be shown using standard arguments, 
it is the eigenfunction of the first excited state of the 
multicomponent Schrodinger-type equation (^7j ) (except 
for a nongcncric case where the components of ip sepa- 
rate) . Since the eigenvalue problem (|27j) is Hermitian for 
B = 0, all eigenvalues A„ with n > 2 are positive, 
the eigenvalue of the ground state is negative, Ao < O.I 

We are not aware of the oscillation theorem for non- 
Hermitian problems. However, since ipi(r) oc rb(r) is 
an eigenfunction with zero eigenvalue even for B =/= 0, as 
B increases from zero, the eigenvalue Ai does not merge 
with other real eigenvalues to form a pair of complex 
conjugate eigenvalues. Therefore, pairs of complex con- 
jugate eigenvalues will be only formed from A n that were 
positive for B = 0. The negative root Ao will remain real 
and negative. In principle, as a result of coalescence of 
complex conjugate eigenvalues, there may emerge pairs 
of negative real eigenvalues. However, the total number 
of negative real eigenvalues will be odd. 



1. Eigenvalues near the potential well 

As an illustration, we consider the eigenvalue problem 
near r we n for the action functional (|I|). Here, Eq. j27| ) 
becomes linear, and the eigenfunctions ■0n( T ) can be 
sought in the form of linear combinations of exp(±iw n T), 
with Lo n = 2nn/(3. The eigenvalues are obtained from 
the equation 



det 



(muj 2 - \ nv )hl + mQ 2 kl 



B, 



0, (30) 



where mfijy = [d 2 U/drkdri] rmU , and B is the magnetic 
field at r we n. The subscript v enumerates the eigenvalues 
A for a given Matsubara frequency. 

If, for example, B is pointing along a principal axes of 
the tensor fijy (say, the axes 1), then we have \ n \jm = 



fif , and 



m 1 A„ 2 ,3 



0J„ 



1 



±- 



1 

f 2 

[(n 



(31) 



where fl 2 , > are the principal values of the tensor 
fij^. Clearly, the eigenvalues A„2,3 are complex conju- 
gate pairs, for large enough uj 2 lu 2 . 

Eq. (31) shows explicitly also how pairs of complex 
eigenvalues emerge with varying magnetic field as a re- 
sult of merging of adjacent real eigenvalues, as discussed 
for the general case in Appendix M. 



B. The prefactor 



We are now in a position to calculate the prefactor in 
the partition function Z. The standard step is to expand 
the deviation Sr(r) of the integration path in (|2^) from 
the extreme trajectory r wc n or i"b in terms of the eigen- 
functions xp n of the corresponding eigenvalue problem, 
<5r(r) = c n'4'n( T )- With account taken of the orthogo- 
nality condition (|29|), the increment SSe of the Euclidean 
action related to the deviation of the trajectory <5r then 
becomes SSe = ^2,^^/2. 

The above expansion assumes that the set {ip n } is 
complete. The completeness is known for B = 0, where 
the eigenvalue problem (^) is Hermitian. As B changes, 
the number of states does not change. From the orthogo- 
nality condition (|2^), none of the wave functions becomes 
a linear combination of other wave functions. This makes 
us believe that the functions rp n form a complete set even 
for B ^ and justifies the above expansion. 

The path integral (^3|) can be obtained as a limit 
N — > oo of integrals over dr(rfc) at discretized instants 
of time Tfc = fcAr, At = j3/N. In the standard way, 
we change to integration over dc n . Because of the or- 
thogonality relation j2S|), the determinant det[i/> n (rfc)] 
of the transformation of variables is real and is equal 
to ±(Ar) Ar ^ 2 . Integration of exp(— SSe) over dc n gives 
const x Y[ n A,T 1/2 - 
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Let us now consider the contribution to the partition 
function Z wc \\ from trajectories close to the potential 
minimum r we ii . The corresponding eigenvalues are 
either positive or belong to complex conjugate pairs, cf. 
(|3l|). Therefore Z wc \\ is real. Since S^r^n] = 0, there 
is no exponentially small factor in -Zwcii- This term gives 
the partition function for low-lying intrawell excitations 
in the presence of the magnetic field. 

In evaluating the contribution from paths close to 
the bounce trajectory, special care has to be taken of the 
eigenvalue A^ = 0. A standard analysifi'El shows that 
integration over dc\ gives the factor (3 in Z\>. The positive 
and complex conjugate eigenvalues A^ give a real posi- 
tive factor in Z\> , whereas the negative eigenvalue A^ (or 
an odd number of negative eigenvalues) make Zb purely 
imaginary. In addition, Z\> contains the exponential fac- 
tor e xpl — 5j;[rb(r)1). Overall, this gives the tunneling 
rate (|2^ ) 

W « 2T|Z b |/Z wcll (x exp{-5 £ [r b (r)]}. (32) 

Eq. (|3^) shows that the instanton technique can be 
used in the presence of a magnetic field in spite of the 
fact that the field breaks time-reversal symmetry. The 
actual calculation is in many respects different from that 
for B = 0. In particular, the bounce trajectory is com- 
plex. The eigenvalues which determine the prefactor 
should be found from a non-Hermitian eigenvalue prob- 
lem. They may be complex, in which case they form 
pairs of complex-conjugate numbers. The bounce trajec- 
tory touches the real escape trajectory at a real point 
rb(/3/2), with a finite real velocity i"b(/3/2). However, 
from our general WKB analysis of tunneling, it follows 
that this is not the point where the particle "shows up" 
as a semiclassical wave packet. 



VII. CONCLUSIONS 

In conclusion, the problem of tunneling in a magnetic 
field can be solved in the semiclassical limit by analyzing 
the Hamiltonian trajectories of the particle in complex 
phase space and time. The boundary conditions are de- 
termined by the intrawell wave function and its analytic 
continuation. This approach allows one to find both the 
tunneling exponent and the tail of the wave function of 
the localized state. It does not require to consider either 
a part of the potential or the magnetic field as a per- 
turbation, and it can be applied to a three-dimensional 
potential of a general form. 

The escape rate is generally exponentially smaller than 
the probability for a particle to reach the boundary of 
the classically accessible range U(r) = E. The escaped 
particle "shows up" from the tunneling barrier with finite 
velocity and behind the surface U(r) = E. The connec- 
tion of the decaying and propagating waves occurs on 
caustics of the set of the Hamiltonian trajectories, where 



the action is branching. The caustics are complex sur- 
faces in 3D space. In the presence of a magnetic field, 
they go through real space along lines (instead of sur- 
faces, for B = 0). 

An interesting feature of tunneling in a magnetic field 
is the occurrence of a switching surface, where there 
merge different WKB branches of the wave function. 
The slope of the logarithm of the wave function sharply 
changes at the switching surface, from its value on one 
of the branches to that on the other branch. The es- 
caped particle first shows up as a propagating semiclas- 
sical wave packet on the switching surface. It happens 
where the classical escape trajectory crosses the switch- 
ing surface. 

The switching surfaces are observable via experiments 
in which the electron density is measured, although such 
experiments arc extremely hard to do. For tunneling 
from 2D electron systems in heterostructures, one can 
think of NMR experiments with samples that contain 
delta-doped layers of nuclear spins (of an isotope that 
differs from that in the host material). If the nuclei are 
sufficiently far from the electron layer, they will detect 
the local electron density on the tail of the wave function 
and its variation with varying fields. One could also use 
light-scattering measurements. 

Switching between branches of the WKB wave function 
for B ^ is similar to the switching between different 
branches of the probability distribution in classical sys- 
tems away from thermal equilibrium. Such systems lack 
time-reversal symmetry, as do also quantum systems in 
a magnetic field. The tail of the classical distribution is 
formed by infrequent fluctuations. Fluctuational paths 
to a given state from the equilibrium position (attrac- 
tor) form a narrow tube centered at the most probable 
path. This path is given by a solution of the variational 
problem of finding the maximum of the logarithm of the 
probability distributior£jEj. In many cases of physical 
interest the corresponding Euler equations are similar to 
Eqs. (^|). However, in contrast to the tunneling problem, 
classical optimal paths can be observedcJ. Switching sur- 
faces in the phase space of fluctuating nonequilibrium 
systems separate areas reached along topologically dif- 
ferent optimal pathscS a ) . They have been seen in analog 
simulationsE3( b ) . 

It follows from the results of this paper that, for poten- 
tial wells which are parabolic near the minimum, even in 
the presence of a magnetic field one can still use the in- 
stanton technique in order to find the escape rate. How- 
ever, the bounce trajectory, which gives the tunneling 
exponent, is now complex. Also in contrast to the B = 
case, the evaluation of the prefactor requires solving a 
non-Hermitian boundary value problem, which generally 
has pairs of complex conjugate eigenvalues. 

The results for the model of tunneling from a strongly 
correlated 2D electron system illustrate the general con- 
clusions about tunneling in a magnetic field. They show 
that the developed method allows us to find the tunnel- 
ing rate and the wave function in a generic system which 
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does not have any special symmetry. They also confirm 
the general conclusions about the structure of singulari- 
ties related to the branching of the WKB wave function, 
and the occurrence of the switching line. We found that 
the tunneling rate in the magnetic field is highly sensi- 
tive to the in-plane electron dynamics and exponentially 
increases when electrons are more strongly confined in 
the plane. It also increases if electrons in the 2D layer 
can adjust to the tunneling electron and thus decrease 
the potential barrier. 

This research was supported in part by the NSF 
through Grant No. PHY-0071059. 



APPENDIX A: THE MANY-ELECTRON 
HAMILTONIAN 

A simple and important model which allows us to con- 
sider the effect of electron correlations on tunneling from 
a 2D electron system is the model of a Wigner crys- 
tal. In this model, the in-planc electron motion is small- 
amplitude vibrations about equilibrium positions. Be- 
cause of strong correlations, exchange effects are not im- 
portant, and the tunneling electron can be identified. Its 
tunneling motion is affected by the interaction with other 
electrons. 

We will assume that the equilibrium in-plane position 
of the tunneling electron is at the origin. Then, in the 
presence of a magnetic field B parallel to the electron 
layer, the full Hamiltonian is of the form 



H=^ + U (z)+H v + H B , 
2m 



(Al) 



with 



H v = - ^2 [ m l PkjP-kj + mo^j-Ukj-u-kj-] (A2) 



and 



H B = -muo 2 c z 2 - LUczN-^^iB x p kj ] z 



k..; 



+U int (z, {u ki }). 



(A3) 



Here, pkj, Uk?, and Wkj are the 2D momentum, dis- 
placement, and frequency of the WC phonon of branch 
j (j = 1,2) with a 2D wave vector k. The in-plane 
momentum of the tunneling electron is N^ 1 / 2 Pkj for 
B = 0. The term Uq(z) describes the tunneling barrier 
[cf. Eq. (^)] for the electron at the origin provided all 
other electrons are at their in-plane lattice sites. 

The term Hb couples the out-of-plane tunneling mo- 
tion to lattice vibrations. The problem of many-electron 
tunneling is thus mapped onto a familiar problptt-of a 
particle coupled to a bath of harmonic oscillatoralSE3. A 



part of the coupling is due purely to the magnetic field. 
Another part comes from the term U m t, which describes 
the change of the tunneling barrier because of electron 
vibrations. Its simplest form is given by the lowest-order 
term of the expansion of the electron energy, 



U in t(z, {U kj }) = 



g-kjU ki 



(A4) 



where gkj are coupling constants [for electrons on a thick 
helium filrrcJ the major term in U- ln t is cx z 2 ]. The cou- 



pling ( A4 ) leads to lowering of the tunneling barrier as a 
result of appropriate displacements of the electrons sur- 
rounding the tunneling electron. 

The major effect on tunneling comes from high- 
frequencii-|in-plane vibrations, which have large density 
of statesEy. Therefore it is not unreasonable to use the 
Einstein model of the Wigner crystal, in which all vibra- 
tions have same freque ncy u ip. Then, except for the term 
Ui n t, the Hamiltonian (Al) becomes a sum of Hamilto- 
nians of confined in the plane noninteracting electrons. 
The Hamiltonian of the tunneling electron has the form 
(^), with the potential U(r) given by <^j. 

For B = 0, the out-of-plane motion of the tunneling 
electron is coupled only to in-plane displacements of other 
electrons. In the Einstein model, it means that the out- 
of-plane motion is decoupled from the in-plane motion 
of the tunneling electron itself. Instead, it is coupled 
to an in-planc oscillator with the coordinate given by a 
(totally symmetric) linear combination of displacements 
of the other electrons. T his m aps the problem onto the 
problem discussed in Sec. VB, with the in-plane electron 
coordinate x in Eq. ( |l§| ) corresponding to the coordinate 
of this oscillator, and with /i being a linear combination 
of the (weighted) coefficients g kj . Of course, for B ^ 
this mapping no longer applies, and extra degrees of free- 
dom have to be taken into account. Yet we expect that 
the model ([l8]) catches much of the qualitative features 
of many-electron tunneling. 



APPENDIX B: EMERGENCE OF COMPLEX 
EIGENVALUES 

In this appendix we consider how, with the varying 
control parameter, two real eigenvalues merge and then 
become complex. Near this transition, the eigenvalues 
can be sought by perturbation theory. We start with a 
value of B = Bo (we can also use another control pa- 
rameter), where the given adjacent eigenvalues A n , A m 
are close to each other and are real. For small |(5B| = 
B — Bo |, the functional T is close to its value for Bo, 
T « in (p7|). To first order in ST, the eigenvalues 

are given by the expressions [A„(B ) + A m (B )]/2 + A±, 
with 
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X ± = -{SF" 



ST 1 



2 



ST" 



ST" 



sx 



AST im ST' 



1/2 



ST nm = & n \ST\iPj. 



(Bl) 



with the wave functions calculated for Bo, and with 
5\ = X m — X n for B = Bq. 

Because of the symmetry (p8[), the matrix elements of 
ST in (Bl) are real. However, the product ST nm 8T mn 
does not have to be positive, and in fact we are interested 
in the case where it is negative. In this case, instead of 
level anticrossing, we have the dependence of the eigen- 
values on the distance SX shown in Fig. 11. In the gap 
the eigenvalues are complex conjugate. 



X 



+ 



NX 
\\ 


// 
/ / 


) > 

//0 


< ( 



FIG. 11. The dependence of the shifted eigenvalues X± on 
the distance 5X between the eigenvalues for B = Bo (schemat- 
ically). We count SX off from ST" n - &T mm , and A± from 
(ST nn + 5T mm )/2. 

We note that the control parameter in the physical sys- 
tem is not SX, and it may be more interesting to look at 
the eigenvalues as functions of the amplitude of ST nm . 
Their behavior is similar to what is shown in Fig. 11, if 
the diagonal and off-diagonal matrix elements depend on 
the control parameter in the same way. Otherwise, once 
the eigenvalues become complex with changing control 
parameter, they do not have to become real again, as is 
the case for the eigenvalues given by Eqs.(||l]) as functions 
of ui c . We note that there is also an opposite process of 
merging of complex conjugate eigenvalues, which is also 
described by Eq. (Bl). 
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